Disappointing model for ultrahigh-energy cosmic rays 
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Abstract. Data of Pierre Auger Observatory show a proton-dominated chemical composition 
of ultrahigh-energy cosmic rays spectrum at (1 — 3) EeV and a steadily heavier composition 
with energy increasing. In order to explain this feature we assume that (1 — 3) EeV protons are 
extragalactic and derive their maximum acceleration energy, E^^^ ~ 4 EeV, compatible with 
both the spectrum and the composition. We also assume the rigidity-dependent acceleration 
mechanism of heavier nuclei, = Z x E^^^ . The proposed model has rather disappointing 

consequences: i) no pion photo-production on CMB photons in extragalactic space and hence ii) 
no high-energy cosmogenic neutrino fluxes; iii) no GZK-cutoff in the spectrum; iv) no correlation 
with nearby sources due to nuclei deflection in the galactic magnetic fields up to highest energies. 



Spectra and chemical compositions of ultrahigh-energy [E >1 EeV) cosmic rays (UHECR) 
measured by two largest detectors, High Resolution Fly's Eye (HiRes) [1] and Pierre Auger 
Observatory (PAO) [2], are significantly different. 

The HiRes data show pure proton composition OH], confirming such signatures of their 
propagation through CMBR as the GZK cutoff [5ll3 and the pair-production dip [71-111]. 
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Figure 1. PAO data [EHU] on Amax(-E') (left panel) and on RMS(Amax) (right panel). Lines 
for protons and Iron are according to QGSJET model |15j . 



The PAO data, on the contrary, strongly favor the nuclei composition getting progressively 
heavier at £" ~ (4 — 40) EeV. This feature, in terms of energy dependence of EAS development 



maximum in atmosphere, X^i^^iE), and r.m.s. of this observable, RMS(Xinax)) is clearly seen in 
Fig. [TJ The data also suggest that the nucleus charge number Z changes smoothly in sources. 

Here we demonstrate that the simple, but disappointing for future experiments, model [16] 
can naturally explain both energy spectrum and mass composition observed by the PAO. 

The basic assumption of the model is the proton composition of UHECR spectrum at 

~ (1 — 3) EeV, the feature supported both by PAO and HiRes. Two more assumptions 
are that these protons are extragalactic and that acceleration of primary nuclei in sources is 
rigidity- dependent, i.e. that E^^^ = Z x Eq, where Eq is a universal energy to be determined 
from data; Z is a nucleus charge number. 

In order to determine the maximum acceleration energy of protons, E^'^^ = Eq, let us 
calculate the extragalactic diffuse proton flux, assuming the power-law generation spectrum 
Qg{E) oc E~'^B with Smax = Eq, and normalize it by the PAO flux at (1 — 3) EeV. Varying 7^ 
in the range 2.0 — 2.8, the maximum value of Eq allowed by the PAO mass composition (see 
Fig. [1]) and energy spectrum (see Fig. [2|) may be obtained. 
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Figure 2. Calculated proton spectra compared to the combined PAO spectrum for different 
E^^^. Extreme cases 7^ = 2.8 and 2.0 are shown in the left and right panels, respectively. 

In our calculations a homogeneous distribution of sources with no cosmological evolution 
(m = 0) was assumed; the highest redshift of sources Zmax = 4. As a criterion of contradiction 
an excess of calculated proton flux at £^ ~ (4 — 5) EeV was chosen. The contradiction has 
different character for different values of 7^. 

For steep source generation functions with 7^ ~ 2.6 — 2.7 the shape and flux of the PAO 
spectrum may be described by E^'^^ 10^*^ — 10^^ eV; the contradiction occurs only in data on 
mass composition. The extreme case, given by 7^ = 2.8, is displayed in the left panel of Fig. [2J 

For flat generation spectra (see the extreme case of jg = 2.0 in the right panel of Fig. [2]) the 
contradiction is very pronounced. For E^^^ = 5 EeV the calculated proton flux exceeds the 
observed one even at -E ~ 2 EeV. 

It is clear that with some redundancy E™^^ ~ (4 — 6) EeV for all 2.0 < 7g < 2.8. 

An influence of possible intergalactic magnetic fields on proton spectrum calculated in a 
diffusive model is shown in the left panel of Fig. [3l Here 7^ = 2.3, which might be the case 
for acceleration by relativistic shocks. The Kolmogorov diffusion in turbulent magnetic field 
with basic scales {Bc,lc) = (1 nG, 1 Mpc) was assumed (see [TTllIS]) and distances between 
sources were d ~ 40 Mpc. The analysis of proton maximum energy of acceleration gives 
again Eq = E^^^ = 4 EeV, in a rough agreement with the analysis made for homogeneous 
distribution of sources. The account for diffusion brings to the flattening of the proton spectrum 
at E < 1 EeV, seen in Fig. [3] as a 'diffusive cutoff', which provides a transition from the steep 
galactic spectrum, most probably composed of Iron, to the flat spectrum of extragalactic protons. 
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Figure 3. Left panel: Comparison of calculated proton spectra with the combined PAO 
spectrum for 7^ = 2.3 and diffusive proton propagation. The cutoff at E^^^ = 4 EeV is needed 
to avoid the contradiction with data at £^ > 3 EeV. Right panel: The energy spectrum in two- 
component model with protons and Iron nuclei with jg = 2.0 and £'max = 4 x Z EeV. The Iron 
nuclei spectrum is calculated for homogeneous distribution of the sources. 



The basic feature of the PAO mass composition, the progressively heavier composition with 
energy increasing, is guaranteed in our model by the rigidity-dependent maximum energy of 
acceleration: at energy higher than Z x E^^^ nuclei with charge Z' < Z disappear, while 
heavier nuclei with larger Z survive. Starting from E^^^ ~ (4 — 6) EeV, the higher energies are 
accessible only for nuclei with progressively larger values of Z. 

Let us now consider a two-component model, with only protons and Iron nuclei being 
produced in sources with generation index 7^ = 2.0 and the maximum acceleration energy 
-E'max = 4: X Z EeV, shown in the right panel of Fig. [3l The primary Iron nuclei spectrum 
is calculated as in |19ll20j for homogeneous distribution of sources. One may notice that the 
calculated spectrum of Iron describes well the cutoff in the PAO spectrum. This steepening is 
caused by the photo-disintegration of Iron nuclei. 

To agree with the mass composition of PAO, the Iron spectrum in Fig. [3] must have a low- 
energy cutoff at E < (20 — 30) EeV. Most naturally it is produced as a 'diffusive cutoff' which 
appears in models with lattice- located sources due to magnetic horizon [21]. Such cutoffs are 
shown in Fig. [Hfor three different sets of parameters Bc,lc,d. The beginning of this cutoff Ec 
for Iron nuclei is Z = 26 times higher than for protons, i.e. Ec ~ 2.6 x 10^^ eV, which has a 
reasonable physical meaning. The gap between 2 EeV and 26 EeV is expected to be filled by 
intermediate nuclei. To provide a smooth RMS(Xinax) curve seen in Fig. [H there are many free 
parameters, e.g. arbitrary fractions of nuclei accelerated in distant sources. 

The predictions of our model are very disappointing for the future detectors. Really, the 
maximum acceleration energy -Emax ~ (100 — 200) EeV for Iron nuclei implies the energy per 
nucleon Ep < E^ax/^ ~ (2 — 4) EeV, well below the GZK cutoff for epochs with z < 15. 
Therefore, practically no cosmogenic neutrinos can be produced in collisions of protons and 
nuclei with CMB photons. Correlation with UHECR sources also is absent due to deflection of 
nuclei in the galactic magnetic fields. The lack of correlation in the model is strengthened by 
the dependence of the maximum energy on Z. 

The signatures of the 'disappointing model' for the PAO detector are the mass-energy relation, 
already seen in the elongation curve Xj^ax{E), and transition from galactic to extragalactic 
cosmic rays below the characteristic energy Ec ^ 1 EeV. 

There are some uncertainties in the model presented above. The most important one relates 
to estimates of E^^^. It is determined by the lowest energy where PAO data become inconsistent 
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Figure 4. As in the right panel of Fig. 3, but 
with the "diffusion cutoff' introduced for three 
different sets of parameters Bc,lc,d. The gap 
between 3 EeV and E'cut (beginning of 'diffusive 
cutoff') is expected to be filled by intermediate 
nuclei with 2 < Z < 25. 



with the proton composition (the 6th low- 
energy bin of the PAO data in Fig. [T|). If 
this energy increases, E'^^^ increases, too. 
The model collapses when the allowed E'!^'^^ 
reaches e.g. (50 - 100) EeV. 

Another case is given by the mass 
composition beinglight nuclei starting right 
from 1 EeV [18j . The cosmological evolution 
of sources are not included in our calculations; 
since this effect slightly decreases E^^^, it 
is not needed to be taken into account. In 
principle, it is also possible that the EeV 
protons detected by PAO are secondary ones, 
i.e. those produced in photo-dissociation of 
primary nuclei in collisions with CMBR and 
extragalactic IR/UV photons. However, in 
fact, as it was demonstrated in [20l[22], the 
flux of secondary protons in the EeV range is 
always smaller than the sum of primary and 
secondary nuclei fluxes. 
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